Gestational age acceleration is associated with epigenetic biomarkers of prenatal physiologic stress exposure.

BACKGROUND: Physiological maternal stress response, such as imbalance in the glucocorticoid pathway and immune system seems to be mediated by DNA methylation (DNAm) and might translate intrauterine stress exposures into phenotypic changes in a sex-specific manner. DNAm in specific sites can also predict newborn gestational age and gestational age acceleration (GAA). GAA occurs when the predicted biological age is higher than the chronological age. In adults, poor health outcomes related to this deviance are well documented and raise questions for the interpretation and prediction in early stages of life. Boys seem to be more vulnerable to intrauterine stress exposure than girls; however, the mechanisms of adaptive sex-specific responses are still unclear. We hypothesize that intrauterine stress exposure is associated with GAA and could be different in boys and girls if inflammatory or glucocorticoid pathways exposure is considered. RESULTS: Using the Western Region Birth Cohort (ROC-São Paulo, Brazil) (n = 83), we calculated DNAm age and GAA from cord blood samples. Two epigenetic risk scores were calculated as an indirect proxy for low-grade inflammation (i-ePGS) and for glucocorticoid exposure (GES). Multivariate linear regression models were applied to investigate associations of GAA with prenatal exposures. The i-ePGS and GES were included in different models with the same co-variates considering sex interactions. The first multivariate model investigating inflammatory exposure (adj. R2 = 0.31, p = < 0.001) showed that GAA was positively associated with i-ePGS (CI, 0.26-113.87, p = 0.049) and negative pregnancy-related feelings (CI, 0.04-0.48 p = 0.019). No sex interaction was observed. The second model investigating glucocorticoid exposure (adj. R2 = 0.32, p = < 0.001) showed that the higher was the GAA was associated with a lower the lower was the GES in girls (CI, 0.04-2.55, p = 0.044). In both models, maternal self-reported mental disorder was negatively associated with GAA. CONCLUSION: Prenatal epigenetic score of exposure to low-grade inflammatory was a predictor of GAA for both sexes. Glucocorticoid epigenetic score seems to be more important to GAA in girls. This study supports the evidence of sex-specificity in stress response, suggesting the glucocorticoid as a possible pathway adopted by girls to accelerate the maturation in an adverse condition.

Sex differences in gene regulatory networks during mid-gestational brain development.

Neurodevelopmental disorders differ considerably between males and females, and fetal brain development is one of the most critical periods to determine risk for these disorders. Transcriptomic studies comparing male and female fetal brain have demonstrated that the highest difference in gene expression occurs in sex chromosomes, but several autossomal genes also demonstrate a slight difference that has not been yet explored. In order to investigate biological pathways underlying fetal brain sex differences, we applied medicine network principles using integrative methods such as co-expression networks (CEMiTool) and regulatory networks (netZoo). The pattern of gene expression from genes in the same pathway tend to reflect biologically relevant phenomena. In this study, network analysis of fetal brain expression reveals regulatory differences between males and females. Integrating two different bioinformatics tools, our results suggest that biological processes such as cell cycle, cell differentiation, energy metabolism and extracellular matrix organization are consistently sex-biased. MSET analysis demonstrates that these differences are relevant to neurodevelopmental disorders, including autism.

Sex differences in gene regulatory networks during mid-gestational brain development

Neurodevelopmental disorders differ considerably between males and females, and fetal brain development is one of the most critical periods to determine risk for these disorders. Transcriptomic studies comparing male and female fetal brain have demonstrated that the highest difference in gene expression occurs in sex chromosomes, but several autossomal genes also demonstrate a slight difference that has not been yet explored. In order to investigate biological pathways underlying fetal brain sex differences, we applied medicine network principles using integrative methods such as co-expression networks (CEMiTool) and regulatory networks (netZoo). The pattern of gene expression from genes in the same pathway tend to reflect biologically relevant phenomena. In this study, network analysis of fetal brain expression reveals regulatory differences between males and females. Integrating two different bioinformatics tools, our results suggest that biological processes such as cell cycle, cell differentiation, energy metabolism and extracellular matrix organization are consistently sex-biased. MSET analysis demonstrates that these differences are relevant to neurodevelopmental disorders, including autism.

Initial findings of striatum tripartite model in OCD brain samples based on transcriptome analysis.

Obsessive-compulsive disorder (OCD) is a psychiatric disorder characterized by obsessions and/or compulsions. Different striatal subregions belonging to the cortico-striato-thalamic circuitry (CSTC) play an important role in the pathophysiology of OCD. The transcriptomes of 3 separate striatal areas (putamen (PT), caudate nucleus (CN) and accumbens nucleus (NAC)) from postmortem brain tissue were compared between 6 OCD and 8 control cases. In addition to network connectivity deregulation, different biological processes are specific to each striatum region according to the tripartite model of the striatum and contribute in various ways to OCD pathophysiology. Specifically, regulation of neurotransmitter levels and presynaptic processes involved in chemical synaptic transmission were shared between NAC and PT. The Gene Ontology terms cellular response to chemical stimulus, response to external stimulus, response to organic substance, regulation of synaptic plasticity, and modulation of synaptic transmission were shared between CN and PT. Most genes harboring common and/or rare variants previously associated with OCD that were differentially expressed or part of a least preserved coexpression module in our study also suggest striatum subregion specificity. At the transcriptional level, our study supports differences in the 3 circuit CSTC model associated with OCD.

Assessment of complementarity of WGCNA and NERI results for identification of modules associated to schizophrenia spectrum disorders.

Psychiatric disorders involve both changes in multiple genes as well different types of variations. As such, gene co-expression networks allowed the comparison of different stages and parts of the brain contributing to an integrated view of genetic variation. Two methods based on co-expression networks presents appealing results: Weighted Gene Correlation Network Analysis (WGCNA) and Network-Medicine Relative Importance (NERI). By selecting two different gene expression databases related to schizophrenia, we evaluated the biological modules selected by both WGCNA and NERI along these databases as well combining both WGCNA and NERI results (WGCNA-NERI). Also we conducted a enrichment analysis for the identification of partial biological function of each result (as well a replication analysis). To appraise the accuracy of whether both algorithms (as well our approach, WGCNA-NERI) were pointing to genes related to schizophrenia and its complex genetic architecture we conducted the MSET analysis, based on a reference gene list of schizophrenia database (SZDB) related to DNA Methylation, Exome, GWAS as well as copy number variation mutation studies. The WGCNA results were more associated with inflammatory pathways and immune system response; NERI obtained genes related with cellular regulation, embryological pathways e cellular growth factors. Only NERI were able to provide a statistical meaningful results to the MSET analysis (for Methylation and de novo mutations data). However, combining WGCNA and NERI provided a much more larger overlap in these two categories and additionally on Transcriptome database. Our study suggests that using both methods in combination is better for establishing a group of modules and pathways related to a complex disease than using each method individually. NERI is available at: https://bitbucket.org/sergionery/neri.

Putative contributions of the sex chromosome proteins SOX3 and SRY to neurodevelopmental disorders.

The male-biased prevalence of certain neurodevelopmental disorders and the sex-biased outcomes associated with stress exposure during gestation have been previously described. Here, we hypothesized that genes distinctively targeted by only one or both homologous proteins highly conserved across therian mammals, SOX3 and SRY, could induce sexual adaptive changes that result in a differential risk for neurodevelopmental disorders. ChIP-seq/chip data showed that SOX3/SRY gene targets were expressed in different brain cell types in mice. We used orthologous human genes in rodent genomes to extend the number of SOX3/SRY set (1,721). These genes were later found to be enriched in five modules of coexpressed genes during the early and mid-gestation periods (FDR < 0.05), independent of sexual hormones. Genes with differential expression (24, p < 0.0001) and methylation (40, p < 0.047) between sexes were overrepresented in this set. Exclusive SOX3 or SRY target genes were more associated with the late gestational and postnatal periods. Using autism as a model sex-biased disorder, the SOX3/SRY set was enriched in autism gene databases (FDR ≤ 0.05), and there were more de novo variations from the male autism spectrum disorder (ASD) samples under the SRY peaks compared to the random peaks (p < 0.024). The comparison of coexpressed networks of SOX3/SRY target genes between male autism and control samples revealed low preservation in gene modules related to stress response (99 genes) and neurogenesis (78 genes). This study provides evidence that while SOX3 is a regulatory mechanism for both sexes, the male-exclusive SRY also plays a role in gene regulation, suggesting a potential mechanism for sex bias in ASD.